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Abstract 

The method of unitary transformations generates five classes of leading-order reducible 
chiral four-nucleon interactions which involve pion-exchanges and a spin-spin contact-term. 
Their first-order contributions to the energy per particle of isospin-symmetric nuclear matter 
and pure neutron matter are evaluated in detail. For most of the closed four-loop diagrams 
the occurring integrals over four Fermi-spheres can be reduced to easily manageable one- or 
two-parameter integrals. One observes substantial cancelations among the different contri¬ 
butions arising from 2-ring and 1-ring diagrams. Altogether, the net attraction generated by 
the chiral four-nucleon interaction does not exceed values of —l.SMeV for densities p < 2pQ. 


1 Introduction and summary 

According to their modern description and construction in chiral effective field theory, nuclear 
forces are organized in a hierarchical way una. For generic few- and many-body observables 
the contributions arising from two-nucleon interactions are larger than those from three-nucleon 
forces, and the latter are again more important than possible corrections due to four-body forces. 
After a renormalization-group evolution to lower resolution scale such chiral low-momentum in¬ 
teractions exhibit desirable convergence properties in perturbative calculations of many-nucleon 
systems and infinite nuclear matter [3llll|5l|6]. In particular, the inclusion of the leading-order chi¬ 
ral three-nucleon interaction (consisting of contact, Ivr-exchange, and 27r-exchange components) 
is essential in order to achieve saturation of isospin-symmetric nuclear matter in this framework. 
The calculation of the sub-leading chiral three-nucleon interaction, built up by many pion-loop 
diagrams etc., has been completed in ref. [7]. By employing chiral two-, three-, and four-nucleon 
interactions up to order N^LO the equation of state of neutron matter at zero temperature has 
been studied in detail by the Darmstadt group [HI 19]. A good convergence pattern has been 
observed for the second and third order particle-particle diagrams with two-nucleon potentials 
only and when including additionally three-nucleon forces in the form of density-dependent two- 
body interactions. Moreover, in ref. [9] the N^LO three-body and four-body contributions to the 
energy per particle of isospin-symmetric nuclear matter have been estimated and large effects 
from the sub-leading chiral three-nucleon force have been found together with small corrections 
from the leading-order chiral four-nucleon interactions. 

An analytical calculation of chiral four-body interactions in nuclear matter and pure neutron 
matter has been performed recently in ref. [10]. The long-range four-nucleon interaction is gen¬ 
erated by pion-exchanges that are combined by the chiral dvr-vertex (off-shell vrvr-interaction) or 
the chiral iViVSvr-vertex. The contributions arising from the 2-ring and 1-ring diagrams related 
to this ’’irreducible” four-nucleon interaction turned out to be very small, with values below 
0.6 MeV for densities p < 0.4fm“^. However, when including the A(1232)-isobar as an explicit 
degree of freedom and counting the delta-nucleon mass splitting A = 293 MeV as a small scale 
(comparable to the pion mass rriT^ or the Fermi momentum kf) a new class of leading-order long- 
range four-nucleon interactions arises. These are mediated by two-fold A(1232)-isobar excitation 
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of nucleons and the exchange of three pions. The analytical calculation of the pertinent 3-ring, 
2-ring, and 1-ring diagrams at four-loop order has lead to a moderately repulsive contribution of 
E{po) = 2.35 MeV at nuclear matter saturation density po = 0.16 fm“^. However, the curve for 
E{p) rises strongly with the density, reaching about 10 times that value at p = 2po. In the case 
of pure neutron matter the same A(1232)-induced four-neutron interactions lead to a repulsive 
contribution that is about half as strong (see Fig. 13 in ref. [T0]h 

The complete set of leading-order four-nucleon interactions in chiral effective held theory has 
been constructed by Epelbaum in refs. mm using the method of unitary transformations. This 
method allows to project the dynamics of the interacting pion-nucleon system into the purely nu¬ 
cleonic subspace relevant for few-nucleon systems at energies below the pion-production thresh¬ 
old. In addition to the ’’irreducible” chiral four-nucleon interaction mentioned above, which 
follows (via the Feynman diagram technique) directly from the vertices of the chiral Lagrangian 
the method of unitary transformations gives rise to eight classes of ’’reducible” four- 
nucleon diagrams. In ref. [12] it is shown that disconnected diagrams always lead to vanishing 
4iV-forces and with the kinematically suppressed pion-energies at a Weinberg-Tomozawa titiNN- 
vertex one obtains that actually three classes (III, VI, and VIII) vanish at leading order. The 
nonvanishing 4V-interactions, grouped into the classes I, II, IV, V, and VII, involve multiple 
pion-exchanges and a spin-spin contact-interaction with coupling constant Ct- The associated 
isoscalar central AW-contact term proportional Cs is absent because it drops out of the commu¬ 
tators arising in the construction of the reducible 4V-interaction. Note that the “irreducible” 
four-nucleon interaction, which constitutes a part of class II, is called -|- V-^ in ref. mi- More¬ 
over, one hnds that in an anti-symmetrized four-neutron state only the (diagonal) matrix-element 
of V“ (class I) is nonvanishing and an explicit expression has been given for it in eq.(Al4) of 
ref. [9]. Analogous expressions for the four-nucleon matrix-elements of the other classes have so 
far not been published. 

The purpose present paper is to evaluate analytically the hve classes (VII, V, IV, II, I) of 
reducible chiral four-nucleon interactions (V*^, V\ V^, V“ in the notation of ref. [TT]) in isospin- 
symmetric nuclear matter and pure neutron matter. In order to facilitate a uniform treatment 
we will hrst introduce in section 2 a general representation of 4V-interactions in terms of a 
product of four (effective) single-nucleon vertices and three ’’propagator”-factors. The latter will 
be ordinary pion-propagators, their squares, or just constants. The hrst-order ”Hartree-Fock” 
contribution to the energy density of nuclear matter is obtained from the connected four-body 
diagram by closing its four nucleon-lines in all possible ways. The actually occurring chiral 
vertices have the convenient property that closing any nucleon-line to itself leads to a vanishing 
spin-trace or isospin-trace. Therefore, one needs to consider here only the closed 2-ring and 
1-ring diagrams generated by the reducible chiral 4V-interactions. These four-loop diagrams are 
then evaluated individually for the hve different classes (in the order VII, V, IV, II, I) in the 
subsections 2.1 up to 2.5. In most cases the occurring integral over the product of four Fermi 
spheres of radius kf can be reduced to an easily manageable one- or two-parameter integral, that 
depends on the dimensionless ratio kf/m.,^. The 1-ring diagrams with crossed pion-exchanges 
related to class II and I allow only for a partial analytical reduction, such that six-parameter 
integrals remain. The numerical results for the energy per particle E{p) as a function of the 
nucleon density p = 2/cj/37r^ are presented and discussed in section 3. We consider two opposite 
values of the spin-spin contact-coupling, Ct = 0.22 fm^ and Ct = —0.45fm^, taken from table I 
in ref. [9] . One observes substantial cancelations between the contributions from different classes. 
In particular, the OT-independent classes II and I provide the largest repulsive and attractive 
contribution, respectively. Altogether, the net attraction generated by the reducible chiral four- 
nucleon interaction stays above values of —1.3 MeV for densities p < 2po = 0.32 fm“^. This 
feature holds for the whole range \Ct\ < 0.5fm^ of the spin-spin contact-coupling Ct- In the 
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Figure 1: Generic form of the ’’reducible” chiral 4N-interaction. The square-box and hlled circle 
symbolize effective single-nucleon vertices. The labels a, b, c denote isospin-indices of exchanged 
pions and i,j denote components of (spin)-vectors. 


case of pure neutron matter, where only class I contributes, the result for the energy per particle 
En{Pn) is determined entirely by a 2 -ring diagram, since the three types of 1 -ring diagrams 
compensate each other almost completely. The corresponding attractive energy per particle goes 
approximately as pH^ and it reaches a value of about —2.1 MeV at pn = 0.4fm“^. 

In summary one hnds that the chiral four-nucleon correlations studied in this work are at 
least one order of magnitude smaller than those provided by the strongly coupled TriVA-system 
[To] with its small mass-gap of 293 MeV. 

2 Diagrammatic calculation of four-body interactions in 
nuclear matter 

In this section we derive analytical expressions for the hrst-order contributions to the energy 
per particle of nuclear matter as they arise from the reducible chiral four-nucleon interactions 
of ref. [12]. In order to facilitate a uniform treatment of the hve different classes (VII, V, IV, 
II, I) we introduce hrst a generic representation of 4A-interactions in terms of a product of 
four (effective) single-nucleon vertices and three ’’propagator”-factors. The latter will be or¬ 
dinary pion-propagators, their squares, or just constants. Fig. 1 shows the intended generic 
form of the reducible chiral four-nucleon interaction, represented as an on-shell scattering pro¬ 
cess AiA 2 A 3 A 4 —)■ The labels a, h, c denote isospin-indices of exchanged pions and 

i,j denote components of (spin)-vectors. These vector-indices are introduced by the spin-spin 
contact-term proportional to Ct, and if present, the left and/or right ”propagator”-factor is a 
constant. The momentum transfers Qa = p{ —pi, Qc = pi—P a and Qb = pi—P 2 + Qa = Ps —pI — Qc 
are given by differences of out- and in-going nucleon momenta. 

From the generic four-body diagram in Fig. 1 the hrst-order ”Hartree-Fock” contribution to 
the energy density pE{p) of nuclear matter is obtained by closing the four nucleon-lines in all 
possible ways. Inspecting the factorized expressions for the hve reducible chiral 4iV-interactions 
written in eqs.(1,5,7,10,19) below, one sees that closing any nucleon-line to itself leads to a 
vanishing spin-trace or isospin-trace. In the case of four neutrons one has either a vanishing spin- 
trace or a zero from the totally anti-symmetric e“^‘’-tensor, since 7 r°-exchange hxed a = b = 3. As 
a consequence of this property, closed 4-ring and 3-ring diagrams vanish trivially and one needs 
to consider here only the closed 2 -ring and 1 -ring diagrams generated by the reducible chiral 
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Figure 2: Closed two-ring diagrams generated by the reducible chiral 4N-interaction. A symmetry 
factor 1/2 occurs, if the square-box and hlled-circle vertex are identical. The diagram at the 
bottom vanishes due to the zero-momentum = 0 of the exchanged central pion. 


4A^-interactions. 

The three possible 2-ring topologies are shown in Fig. 2. One notices that the diagram at 
the bottom vanishes due to the zero-momentum iff, = 0 of the exchanged central pion. At least 
one of the adjoined vertices involves qb linearly. Furthermore, the six possible 1-ring topologies 
are shown in Fig. 3. According to the shape of the dashed (pion) line the diagrams in the first, 
second, and third row are referred to as U-type, Z-type, and X-type, respectively. Note also that 
if the square-box and hlled-circle vertex represent identical vertices (this happens for class VII 
and I) a symmetry factor 1/2 belongs to the 2-ring diagrams and to the Z-type 1-ring diagrams, 
and obviously identical copies of the U-type and X-type 1-ring diagrams need to be discarded. 
After these preparations we can now turn to the evaluation of the individual 1-ring and 2-ring 
diagrams generated by the hve classes of reducible chiral 4V-interactions. 
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Figure 3: Closed one-ring diagrams generated by the reducible chiral 4N-interaction. The dia¬ 
grams in the hrst, second, and third row are referred to as U-type, Z-type, and X-type, respec¬ 
tively. A symmetry factor 1/2 belongs to the Z-type diagrams, if the square-box and hlled-circle 
represent identical vertices. 


2.1 Class VII 


We start (in reverse order) with class VII, which is quadratic in the spin-spin contact-coupling 
Ct- The corresponding factorized expression for this chiral 4V-interaction reads: 
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ch\ 


a\ (^2 X QbY T. 
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(as X qbY Tg 
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( 1 ) 


Note that there is an extra factor 2 in comparison to eq.(3.60) in ref. [12] due to the interchange 
symmetry (12) y-)- (43) between pairs of nucleons. A straightforward evaluation of the left 2-ring 
diagram in Fig. 2 gives the following contribution from class VII to the energy per particle of 
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isospin-symmetric nuclear matter: 


E{p) = 


[ dy x{u^ - y'^){u - xy{2u + x){ - 8xy + 6x 
[271 u) Jo Jo 


X [arctan(2a; -|- 2y) — arctan(2a; — 2y)\ -|- {2x^ — 2y‘^ — 1) In 


1 -I- A[x + y)'^ 
1 -I- 4:{x — yY 


, (2) 


with the dimensionless variable u = kf/rriT^. The nucleon density p = 2fcj/37r^ is related to the 
Fermi momentum kf in the usual way. In order to arrive at a double-integral representation we 
have used the reduction formula: 


(, l*\ I'LL 

dp I dzp^F{pz \/v? — p^{l — z"^)) = I dy {u^ — yYF{2y), 


( 3 ) 


to eliminate a complicated angular integral, and a master formula for integrating F’dpi—P 2 I) over 
two Fermi spheres |pi, 2 | < kf. The latter introduces the weighting-function x‘^{u — xY{2u + x). 
Actually, there exists an analytical solution of the double-integral in eq.(2), namely: 
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Note that in the chiral limit 0 only a term M^(1281n2 — 179)/10 = —9.0277 m^ remains. 

Let us also comment on the features of the other diagrams for class VII: The right 2-ring diagram 
in Fig 2 leads to a vanishing isospin-trace in symmetric nuclear matter, whereas in pure neutron 
matter the left and right 2-ring diagram cancel each other. The U-type and Z-type 1-ring 
diagrams in Fig. 3 cancel each other after taking the spin-trace, and the X-type 1-ring diagram 
leads to a vanishing spin-trace by itself. 


2.2 Class V 

Next we come to class V, which is induced by the Weinberg-Tomozawa TrvriViV-vertex in combi¬ 
nation with the spin-spin contact-coupling Ct- The corresponding factorized expression for this 
chiral 4iV-interaction reads [T2] : 
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and it obviously vanishes for four neutrons, where only 7r®-exchange is possible = 0). The 
evaluation of the X-type 1-ring diagram in Fig. 3 gives the following contribution from class V to 


6 
























the energy per particle of isospin-symmetric nuclear matter: 


E{p) = 


?,CTg\ml 
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Note that the integral j^dy ... could be solved in terms of elementary functions (In, arctan), but 
we refrain from presenting this rather lengthy expression. The same applies to the subsequent 
results given in eqs.(8,9,ll,20,21). Concerning the other diagrams for class V, one hnds that both 
2-ring diagrams in Fig. 2 lead to a vanishing spin- and isospin-trace. Moreover, the U-type and 
Z-type 1-ring diagrams in Fig. 3 cancel each other after taking the spin-trace. 


2.3 Class IV 


Next we come to class VI, which is induced by the spin-spin contact-coupling Ct in combination 
with a ’’contraction” of two ordinary 7riV-vertices. The corresponding factorized expression for 
this chiral 4V-interaction reads 
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Whereas the left 2-ring diagram has a vanishing isospin-trace, one obtains from the right 2-ring 
diagram the following contribution to the energy per particle of isospin-symmetric nuclear matter: 
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By rewriting the spin- and isospin-factors as a product of two anti-commutators one hnds that the 
U-type and Z-type 1-ring diagrams add to zero. On the other hand, the X-type 1-ring diagrams 
in Fig. 3 give a contribution to E{p) of the form: 
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In the case of pure neutron matter the left and right 2-ring diagram cancel each other, likewise do 
this the U-type and Z-type 1-ring diagrams, and the X-type 1-ring diagrams vanish by themselves 
(because for four neutrons the e“^'^-term in is absent). 


2.4 Class II 

Next we come to class II, which is induced by the Weinberg-Tomozawa vrTriViV-vertex in combina¬ 
tion with the ’’contraction”of two ordinary vriV-vertices. The corresponding factorized expression 
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for this 37r-exchange 4iV-interaction reads 
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ml + ^ 

and it obviously vanishes for four neutrons, since = 0. The evaluation of the left 2-ring 
diagram in Fig. 2 gives the following contribution to the energy per particle of isospin-symmetric 
nuclear matter: 
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whereas the right 2-ring diagram produces a vanishing spin-trace. In the case of the U-type 
1-ring diagrams in Fig. 3, three (out of four) Fermi sphere integrals over the pion-propagators 
and momentum-dependent interactions factorize by making use of tensor-contractions. For this 
reason one can represent the corresponding contribution to E{p) as a one-parameter integral of 
the form: 
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defined by vectorial and tensorial Fermi sphere integrals over a pion-propagator (see also eq.(A.l) 
in ref. pro]). In the case of the Z-type 1-ring diagrams in Fig. 3, the two Fermi sphere integrals over 
the pion-propagators associated to the short pion-lines factorize via tensors. For the remaining 
two Fermi sphere integrals the angular part can be carried out, and in this procedure the third 
pion-propagator introduces the logarithmic function: 

l + (x + i /)2 
1 4- (x — yY 

Putting all pieces together, one obtains the following double-integral representation for the con¬ 
tribution of the Z-type 1-ring diagrams from class II to the energy per particle of nuclear matter: 
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In the case of the X-type 1-ring diagrams in Fig. 3 the assignment of momenta to the virtual pions 
is now such that the previous factorization approach does not work anymore. Only the Fermi 
sphere integral associated to the nucleon-line between the square-box and hlled-circle vertex can 
be performed analytically. Hence, the contribution of the X-type 1-ring diagrams from class II 
to the energy per particle of nuclear matter result in the form: 
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with qi = Pi — Ps, q 2 = P 2 — P 3 and ff = Pi + P 2 — Ps, where all momentum-vectors have been 
divided by the pion mass Note that due to rotational invariance it is sufficient to parametrize 
the nine-dimensional integral in eq.(18) by three radii (0 < pi, 2,3 < u), two directional cosines 
(—1 < zi^2 < 1), and one azimuthal angle (0 < (p < 27r). 


2.5 Class I 


Finally we come to class I, which is induced by a pair of ’’contractions” of two ordinary ttN- 
vertices. The corresponding factorized expression for this Svr-exchange 4A^-interaction reads: 
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Note that there is again an extra factor 2 in comparison to eq.(3.33) in ref. [12] due to the 
interchange symmetry (12) t-)- (43) between pairs of nucleons. The evaluation of the left 2-ring 
diagram in Fig. 2 gives the following contributions to the energy per particle of isospin-symmetric 
nuclear matter: 
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and to the energy per particle of pure neutron matter: 
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The neutron density pn is related to the neutron Fermi momentum kn by pn = kl/STr"^. In eq.(21) 
and all following formulas for En{pn) the dimensionless variable u has the meaning u = kn/m-,^. 
Next one notices that for class I the right 2-ring diagram in Fig. 2 produces vanishing spin-trace. 
Furthermore, one obtains from the U-type 1-ring diagram in Fig. 3 the following contributions: 


E{p) = 


9aK 


A{A7rUyu^ /o 


dx 


8GUx) + GUx) Ks{x) + 8Gs{x) + Gt{x) Gt{x)Kt{x) \ , (22) 


9 





















aWn 


En{Pn) = I2(4y^6^3 { [^GI{x) + GUx)]k s{x) + [Gt{x) - AG s{x)]Gt{x) Kt{x)^ , (23) 

where we have introduced two new auxiliary functions: 

/ \ o / \ o / \ 2 + 1 + (m + X)'^ 

Es[x) = 2u — 6 arctanfM + x) — 3 arctanfw — x) -\ -In--- — 

2x 1 + (m — xy 


+ 2m' “ 2) - 3(1 + In 

dehned by a tensorial Fermi sphere integral over a squared pion-propagator. 

On the other hand the Z-type 1-ring diagrams give the following contributions: 
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arising from an angular integral f^^dz... over a squared pion-propagator. 
Finally, the X-type 1-ring diagram leads to contributions of the form: 


E(p) = 


3s‘ml 


d^Pi d^p2 d?po 


{ATrUyU^ J (27r)3 (1 + g-;2)(l + ^-|) 

\pj\<U 


X 3gf g| 




2KT{\f]\) +{qi-q2y 


?2)'-7((fi-g2)(h-gi)(h-92) h 


EniPn) 


alml 

(47r/^)S3 


d^Pi d^p2 d^po 1 

(i+^Ki+W) 


Xi 3(gi ■ ga)' 


+ i7T 


^2 ^2 
Q2 


q 2 f - 3(gi ■ q 2 ) {ff- qi) {ff- ^ 2 ) 


10 




Figure 4: Reducible chiral 4iV-interactions in isospin-symmetric nuclear matter. The dashed 
lines show the individual contributions from the hve classes for Ct = 0.22 fm^, and the full line 
gives their total sum. 

with (fi = Pi - p 3 , q 2 =p 2 - p3 and ff = pi + p 2 - ps- 

This completes our presentation of semi-analytical results for the energy per particle of isospin- 
symmetric nuclear matter E{p) and pure neutron matter E{p) as derived from the reducible chiral 
4A^-interactions. The given formulas allow for an easy and accurate numerical evaluation [1^ of 
these leading-order chiral four-body correlations. 

3 Results and discussion 

We are now in the position to present numerical results. The physical parameters related to 
pion-exchange are: qa = 1-29 (nucleon axial-vector coupling constant), /tt = 92.4MeV (pion 
decay constant), and = 135 MeV (neutral pion mass). For the spin-spin contact-coupling Ct 
we choose two opposite values, Ct = 0.22 fm^ and Ct = —0.45fm^, taken from table I in ref. [9]. 

The contributions to the energy per particle E{p) of isospin-symmetric nuclear matter are 
presented in Fig. 4 for the small CT-value, Ct = 0.22 fm^, in the density region po/4 < P < 2po; 
with po = 0.16 fm“^ the empirical saturation density. The dashed lines (with appropriate labels) 
show separately the hve classes (VII, V, IV, II, I) and the full line gives their total sum. One 
observes that class VII (proportional to C^) is smallest in magnitude and that other classes (V-I-IV 
and II-I-I) cancel each other to a large extent. The largest repulsive and attractive contributions 
are provided by the (Fr-independent classes II and I, respectively. As a consequence of these 
compensations the net attraction stays above —1.3 MeV for all densities p < 0.36 fm“^. The 
analogous results for the larger negative OT-value, Ct = —0.45fm^, are shown in Fig. 5. One 
notices that class VII (~ 0|>) is still smallest in magnitude and that the classes V and IV have 
changed their position in comparison to the sequence in Fig. 4. The net attraction reaches now 
a value of —1.7 MeV at p = 0.36 fm“^. Exploring the whole range \Ct\ < 0.5 fm^ of the spin-spin 
contact-coupling Ct, one hnds that the 4V-attraction in nuclear matter does not exceed values 
of —1.3 MeV for densities p < 2po. 


11 






Figure 5: Reducible chiral 4iV-interactions in isospin-symmetric nuclear matter. The dashed 
lines show the individual contributions from the hve classes for Ct = —0.45fm^, and the full line 
gives their total sum. 



Figure 6: Reducible chiral 4iV-interactions in pure neutron matter. The conributions from 2-ring 
and 1-ring diagrams are shown separately. The U-type, Z-type, and X-type diagrams cancel each 
other almost completely. 

In the case of pure neutron matter only class I (see subsection 2.5) contributes and the 
result for the energy per particle En{pn) comes out parameterfree. Fig. 6 shows by the dashed 
lines the contributions from the 2-ring diagram and the sum of the three 1-ring diagrams in the 
density region 0.04 fm“^ < pn < 0.4 fm“^. One recognizes that the latter part is extremely small. 
This feature has its origin in the almost complete compensation between the contributions from 
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the U-type, Z-type and X-type 1-ring diagrams, which are shown by the three dashed-dotted 
lines in Fig. 6. Interestingly, the net resnlt for En{pn) is thus entirely determined by the 2-ring 
part written in eq.(21). The corresponding attractive energy per particle goes approximately 
as and it reaches the value of about —2.1MeV at the (relatively high) neutron density of 
pn = 0.4 fm“^. 

In passing we note that our numerical results for the chiral four-body correlations in nuclear 
and neutron matter agree perfectly na with those of the Darmstadt group [9] if the (unessential) 
regulator-function is set to /r = 1. 

In summary one can conclude that the chiral four-nucleon correlations studied in this work are 
at least one order of magnitude smaller than those provided by the strongly coupled TriVA-system 
im with its small mass-gap of 293 MeV ~ 1.1kfQ. 
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